Phase transitions in full counting statistics for periodic pumping 
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We discuss the problem of full counting statistics for periodic pumping. The probability generating 
function is usually denned on a circle of the "physical" values of the counting parameter, with its 
periodicity corresponding to charge quantization. The extensive part of the generating function 
can either be an analytic function on this circle or have singularities. These two cases may be 
interpreted as different thermodynamic phases in time domain. We discuss several examples of 
phase transitions between these phases for classical and quantum systems. Finally, we prove a 
criterion for the "analytic" phase in the problem of a quantum pump for noninteracting fermions. 



I. INTRODUCTION 

The problem of full counting statistics (FCS)^is often 
considered in a setup periodic in timer In such a formula- 
tion, the system depends on external parameters varying 
periodically in time, and one is interested in counting 
certain quantized events (typically, the transfer of par- 
ticles between the leads of a contact). Generally, the 
probabilities of different outcomes P n are labeled by in- 
teger indices n and can be combined into the probability 
generating function^ 



+ OG 



X (A)= J2 p ^ 1 



(1) 



For a superposition of statistically independent processes, 
the generating function is given by the product of those 
for each process. Therefore, for a periodic process ex- 
tended in time, this generating function is exponentially 
extensive in time, provided the correlations in time de- 
cay sufficiently rapidly. One therefore usually defines the 
generating function "per period'™ (or "extensive" FCS) 



Xo(A) =cxp 



m^ 1 ^ ln x( A ) 

iV— >-00 1\ 



(2) 



where the FCS xW i s collected over N periods. 

A vast literature is devoted to FCS both in the general 
setup and for periodic processes.— Probably, the most 
interesting class of FCS problems are those of quantum 
charge transfer)^ but there are also discussions of FCS 
in classical stochastic processes (see, e.g., Ref. [S]) and of 
the relation between classical and quantum effects in FCS 
(see, e.g., Sec. 5 of Ref. 4 and Ref. [5]). 

From this immense body of results, one notices that 
FCS in time-extensive problems can be conveniently clas- 
sified by the analytic properties of the FCS generating 
function per period Q. Namely, two main phases can 
be identified: lnxo(^) may either be analytic on the unit 
circle of real A or have a singularity at certain values of 
A. Correspondingly we distinguish two phases of FCS: 
analytic and nonanalytic. 



In the present paper we show that the existence of 
these two phases is a very general feature of FCS: phase 
transitions between them occur both in classical stochas- 
tic models and in quantum systems. We illustrate our 
discussion with several examples: the classical weather 
model and quantum systems of noninteracting fermions. 
For quantum examples, we use the results of our recent 
wor yZm |; show the stability of the analytic phase for a 
wide class of non-interacting fermionic systems, even at 
finite temperature. Finally, we discuss a possible iden- 
tification of the two phases in terms of cumulants and 
conjecture the general form of the asymptotic long-time 
behavior of the FCS in the nonanalytic phase. 



II. TWO PHASES OF FCS IN TIME-PERIODIC 

SYSTEMS 

The generating function ([I]) is periodic, x(A) = x(^ + 
27r), and its Fourier components P n are non- negative. 
Furthermore, the probabilities P n always obey the nor- 
malization condition y^„ P n — 1, and therefore for non- 



negative P n the series (11 
1. Typica 



circle \e lX 



) converges uniformly on the unit 
ly, P n decay sufficiently rapidly 
as functions of n, and then xM is analytic for real values 
of A. 

Out of the above three properties (periodicity, non- 
negativity of probabilities and analyticity) only the first 
one (periodicity) necessarily holds for the generating 
function per period (pi) from its definition. As we shall see 
below, both the non- negativity of the Fourier components 
( "quasiprobabilities" ) and the analyticity may be broken 
for Xo(A). Note that non- negativity of the quasiprobabil- 
ities is not always related to analyticity (e.g., Xo(-M may 
be analytic and still have some negative quasiprobabili- 
ties). 

As we shall see below, most common nonanalityc fea- 
tures of XoW are discontinuities and kink points. In 
those cases, the quasiprobabilities decay algebraically in 
n and oscillate in sign (if the singularity is located at 
A 7^ 0). On the other hand, an analytic Xo(A) corre- 
sponds to an exponential decay of the quasiprobabili- 





FIG. 1: (a) Weather model. Transitions to sunny (rainy) 
weather S (R) are shown by arrows with corresponding prob- 
abilities, (b) Phase diagram of the weather model in the 
(q a , q r ) parameter space. The shaded area corresponds to the 
nonanalytic phase (NA). 



value remains leading for all values of A (analytic phase) 
or the leading eigenvalue switches at some value of A 
(nonanalytic phase). From the explicit formula for the 
eigenvalues 

q r +q s e lX \( q r + q s e lX \ 2 , . ,, 

wi,2 = 2 ± Vl 2 J +( 1 -^-«'-) e 

one can deduce the phase diagram in the (q Sl q r ) coordi- 
nates, (Fig. lb) and find that the nonanalyticity occurs 
at A = it. 



ties (not necessarily non-negative). Thus we may view 
the difference between the analytic and nonanalytic be- 
haviors of Xo(A) as a distinction between two quantum 
phases. As usual for phase transitions, the transition 
to the nonanalytic phase appears only in the thermody- 
namic limit (when one considers the extensive part of the 
generating function). 

We may also remark in passing that there is nothing 
surprising in negative quasiprobabilities in the extensive 
generating function Xo(A). Indeed, this generating func- 
tion appears in our attempt to factorize the full FCS into 
N independent processes |2]). Generally, there is no rea- 
son to expect that such decomposition is possible with 
physical non-negative probabilities, since the system has 
certain correlations between different periods. 



III. CLASSICAL EXAMPLE: WEATHER 
MODEL 

The simplest stochastic process that illustrates the 
analytic-nonanalytic phase transition is the so called 
"weather model" . 9 Suppose that the weather on each day 
can be of two types: rainy or sunny, and that the weather 
on a given day is chosen randomly with the probabilities 
depending on the weather on the previous day. Let the 
probability of a sunny day after a sunny day be q s and 
the probability of a rainy day followed by a rainy day 
be q r (Fig. la). This defines a stochastic process, in 
which the FCS for the number of sunny days can be eas- 
ily calculated.^ The calculation involves two generating 
functions Xs(A) and Xr-(A), which correspond to the con- 
ditions that the last day is sunny or rainy, respectively. 
The evolution of these generating functions is described 
by a linear operator acting in the two-dimensional space. 
The full generating function is then given by 

X (X) = xs(\) + XrW = Ci[wi(A)] w + C 2 [u 2 (\)] N , (3) 

where N is the number of observation days, wi(A) and 
W2(A) are the two eigenvalues of the evolution operator, 
and C\ and C 2 are some coefficients. Therefore 



Xo(A) 



:( Wl (A),u; 2 (A)) 



(4) 



IV. QUANTUM EXAMPLES: 
NON-INTERACTING FERMIONS 

As quantum examples illustrating the analytic- 
nonanalytic transition in time-periodic FCS, we consider 
systems of non-interacting fermions. This class of sys- 
tems has been studied extensively,® with the main result 
given by the so-called Levitov-Lesovik determinant for- 
mula for the characteristic function (TH)) 1 ! 3 ! 11 !!!]]^ ^g p e _ 
riodic case, the same determinant formula can be used to 
find the extensive part of the FCS (pi) by simply imposing 
periodic boundary conditions in time. 3 

For simplicity, we consider here a one-channel quantum 
contact of non-interacting spinless fermions. The trans- 
parency of the contact g is taken to be time-independent, 
but we assume a time-dependent voltage V(t) applied to 
the contact and an arbitrary temperature T. This sys- 
tem was considered in previous works p^l and we can use 
their results to study the analytic-nonanalytic transition 
in a periodic setup. The case of a more general time- 
dependent scattering matrix-^ will be treated in the next 
section. 

The FCS for non-interacting fermions can be expressed 
in terms of the distribution function of effective trans- 
parencies n(p) for single-particle processes.® This func- 
tion determines the jump of the derivative of lnxo(A) 
on the negative real axis of the variable e lX [so that 
V 6(0,1)]: 



Kp) 



2ni 



dplnx(A) 



p—iO 

p+?0 



oiX 



(6) 



where the eigenvalue with the maximal absolute value is 
chosen. Then two phases are possible: either one eigen- 



Since, in the non- interacting case, singularities of Xo(A) 
are allowed only on the negative real axis of e lX (i.e., 
if we take A to be real, at A = tj-) EEl the analytic and 
nonanalytic phases correspond to /i(p = 1/2) = and 
(j,(p = 1/2) y^ 0, respectively. 

From the previous studies of the one-channel contact 
with a time-independent transparency g, we can sketch 
the general behavior of the distribution function /i(p) in 
several special cases (Fig. 2): 

(a) T = 0, V(t) = const. In this case n(p) — aS(p — g) 
or /j,(p) = aS(p — [1 — g]), depending on the polarity of 
V (Fig. 2a). The weight a = Vt/(2tt), where r is the 
conventional period. This case obviously belongs to the 
analytic phase, unless g = 1/2. 



t 



iT=0, V=const 



AA_ 



7>0, V=0 



g or (1-g) 



T ► 

1 (a) 



l-J 



P 

-► 



(l±VT^)/2 



1 (c) 



AA 


7=0, any VffJ 


1 


, 1 — us. 



kJL 



T>0, V=const 



g and (1-g) 



1 (b) 



(1+VT^)/2 (d) 



FIG. 2: Schemes of distributions of effective transparencies 
fi(p) for different examples of FCS in a quantum contact. 
Shaded areas represent a continuous spectrum u(p). Peaks 
at gap edges schematically represent delta-function contribu- 
tions. 



(b) T = 0, arbitrary V(t). In this case, /x(p) is given 
by a superposition of the delta function from the pre- 
vious example and a generally continuous spectrum for 
\p — 1/2] > \g — 1/2 1 (gapped around p — 1/2, except 
for g — 1/2), see Fig. 2b. This form of the spectrum 
(delta function and a continuous part) follows from the 
arguments in Refs. I7IT5|) . In addition, it is proved there 
that the continuous part is symmetric with respect to 
p = 1/2. The existence of the gap around p — 1/2 also 
follows from the theorem proved in the next section. Ex- 
cept for g = 1/2, this case also belongs to the analytic 
phase. 

(c) T > 0, V(t) = 0. In this case, the spectrum /i(p) 
is symmetric, continuous, and gapped around p = 1/2 
(Fig. 2c). The value of the gap can be obtained by a 
direct calculation^ 21 ^ \ p - 1/2| > y/T^g/2 (also pre- 
dicted by the theorem proven in the next section). This 
example belongs to the analytic phase as well. 

(d) T > 0, V(i) = const (without loss of generality, 
we assume V > 0, otherwise we may simply reflect p t— > 
1—p). In this case, a calculation shows that the spectrum 
is continuous and consists of two regions: < p < g and 
(1 + y/T^g)/2 < p < 1 (Fig. 2d). This case belongs 
to the analytic phase, if g < 1/2 and to the nonanalytic 
phase, if g > 1/2. Note that the two regions of support 
of /x(p) overlap if g > 3/4. 

(e) T > 0, arbitrary V(t). In this case, the spectrum 
fi(p) is generally continuous and non-symmetric. As we 
prove in the theorem in the next section, the analytic 
phase is realized for g < 1/2 (which results in the gap 
\p — 1/2 1 > 1/2 — g). The nonanalytic phase is generic 
for g > 1/2. 



V. ANALYTIC PHASE FOR 
NON-INTERACTING QUANTUM SYSTEMS 

Let us now consider a slightly more general setup: a 
one-channel two-lead contact which can be described by 
an instantaneous time-dependent 2x2 scattering matrix 
S(t). We assume a finite temperature T. The character- 




FIG. 3: Spherical cap of height 2go covering the trajectory 
N(i) and the north pole e z . 



istic function xW m this case is given by^l 



X (A) = det 1 + (e lA - l)X 



-iXQ 



(7) 



Here the determinant is taken in the single-particle 
Hilbert space, Q = (1 + a z )/2 is the charge operator in 
one of the leads, and the Hermitian operator X is defined 
as {uf is the Fermi occupation number) 



X = (1 - n F )Q + (npy^S^QSinF) 1 / 2 . 



(8) 



The spectrum of the operator X is real and confined to 
the interval [0,1]. Its eigenvalues can be interpreted as 
effective transparencies of the contact between leads, and 
its spectral density /i(p) = tr S(p — X) determines the 
characteristic function as 

ln X (X) = -iXTrQ+ [ d P fi{p)\n[l + (e tX -l)p] . (9) 



The same calculation applies to the extensive part of the 
generating function xo PO > H one imposes periodic bound- 
ary conditions in timeP 

As mentioned in the previous section, the analytic 
or nonanalytic phase in the periodic setup depends on 
whether fi(p — 1/2) is zero or non-zero. The proof 
of the analytic phase thus amounts to demonstrating 
a gapped region around p = 1/2 for a certain class of 
time-dependencies S(t). Namely, we prove below that 
if the trajectory of S(t) is confined within a certain re- 
gion (spherical cap), then there is a gap in the spectrum 
n[jp) around p = 1/2. Consider the matrix N — S^a z S 
or, equivalently, the vector N — (l/2)Tr<7JV (it is a time- 
dependent vector on the unit sphere^!) . We suppose that 
this vector, at all times, together with the north pole 
N = e z , fit in a certain spherical cap of height 2go < 1 
(see Fig. [3]). Then we can prove that (i(p) — within the 
window \p — 1/2 1 < 1/2 — go. 

Indeed, the above assumption may be expressed math- 
ematically as the existence of a constant unit vector M 
such that 

MN(t) > 1 - 2.9o for all t and Me z >l-2g . (10) 

We further denote M = Ma and assume that np is a 
scalar in the lead space (the same temperature in both 



leads), in particular [tif,M] 
a simple algebra, 



\nF,<J 



0. Then, after 



2X - 1 > - 



-(aM) 2 + { 



2X - l,aM 



} 



= —a + a 



+ a ({M, <j z }{\ - n F ) + nJ 2 {M, N}nJ 2 ^j 



(11) 

This inequality holds for any real a. In particular, we 
can choose a = 1 — 2go > 0. If we now use the fact that 
the eigenvalues of rip are all located between and 1, we 
arrive at the desired inequality 



(^-l/2) 2 >(l/2-. 9o )' 



(12) 



This result implies a gap in (i(p) of the size 1 — 2<?o cen- 
tered around p — 1/2, i.e., the analytic phase. 

The class of trajectories for which our theorem applies 
is quite wide. One particular case are those trajectories, 
for which the instantaneous transparencies g{t) = (1 — 
N(t)e z )/2 are bounded from above by a certain maximal 
transparency go < 1/2. In this case, one may safely 
choose M = a z and prove the analytic phase. 

Interestingly, in the example (c) of the previous section 
(T > and time-independent S 1 ), it is more advantageous 
to tilt the vector M so that it bisects the angle between 
e z and N: this expl ains th e gap yl — g in the spectrum 
/j,(p) in this example! 12 * 14 ! 

Another interesting particular case is T = 0. In this 
case, /i{p) is invariant with respect to global rotations 
of the trajectory N(t) (see Ref. 15). This implies that 
the condition Me z > 1 — 2go does not need to be taken 
into account: it suffices to cover only the trajectory N(t) 
(and not the north pole) by any spherical cap smaller 
than hemisphere, in order to prove the analytic phase. 
In particular, the transparencies g and 1 — g produce the 
same gaps in examples 1 and 2 of the previous section. 

Note that the above theorem may also be extended 
to the case of different temperatures of both leads. In 
this case, we are restricted to M = a z (the condition 
[np,a z ] = is still assumed as the absence of initial 
entanglement of the leads^). 

Finally, the theorem also applies to the multichannel 
prob lem, with the only modification that the conditions 
( fl0| should be replaced by (1/2) {M, N(t)} > l-2g and 
(1/2){M, a z } > 1 — 2g Q , in terms of all eigenvalues. 



VI. DISCUSSION OF RESULTS 

We have identified two phases of FCS, analytic and 
nonanalytic, in terms of the analyticity of the generat- 
ing function per period (pi) at real values of the counting 
parameter A. Although singularities of Xo(A) in all exam- 
ples considered in the present paper appeared at A = n, 
they may generally occur at any values of A. We sug- 
gest that the appearance of those singularities may be 



the most common scenario of the analytic-nonanalytic 
phase transition and ask the natural question: what are 
the fingerprints of such a transition? 

The answer to this question may be deduced from 
analyzing classical models (e.g., the weather model de- 
scribed above). Let us define the staggered average: 
{A{n))^ = ((-l) n A{n))/((-l) n ) for any function of the 
counted events A(n)y^Then the staggered cumulant 



((n 2 )), = (n 2 ) n - (n)l = Hd A ) 2 ln X (A)| A= 



(13) 



has different asymptotics as a function of the number 
of periods N in the two phases. In the analytic phase, 
((n 2 )),r oc N, while in the nonanalytic phase, two eigen- 
values Wi,2(A) contribute, and one finds ((n 2 ))^ ex N 2 . 

In quantum noninteracting systems, the same distinc- 
tion between the two phases can me made if the asymp- 
totic behavior of the generating function x(A) over N 
periods is given in the nonanalytic phase by 



X(A) = C\[ X o(A - 0)] JV + C 2 [ X o(A + 0)] 



iN 



(14) 



where A ± refers to taking the analytic branches be- 
low and above the cut at A = ir. The coefficients C1.2 
are non-extensive (sub-exponential in N) and generally 
depend on A. In some situations (e.g., for transmis- 
sion of non-interacting fcrmions through an open channel 
or for counting one-dimensional fermions on a line seg- 
ment at zero temperature), these coefficients are known 
to have an algebraic time dependence, C1.2 oc 7V 7 ( A=F °), 
and thus contribute logarithmically to cumulants. 1 - In 
fact, in these simple examples, the expression (14) can 



be rigorously proven using the extended Fishcr-Hartwig 
conjecture (theorem)) 17 1 18 1 as the FCS in those cases is 
given exactly by a Toeplitz determinant. 



We may further conjecture that the expression (14 1, 
with Ci,2 algebraically depending on N, should remain 
valid for a larger class of models: for example, for all non- 
interacting fermionic systems p where the FCS problem 
can be reduced to the determinant of a block- Toeplitz 
matrix. It may even be possible that this result is appli- 
cable for interacting systems. These interesting problems 
remain for future study. 

Note that a phase transition in FCS in a particular in- 
teracting system has been recently discussed in Ref. [13 
with the difference that their phase transition reveals 
itself in the ordinary (non-staggered) cumulants and is 
likely connected to a nonanalyticity in Xo(A) at A = (cf. 
A = 7T in our examples). One can also invent other ex- 
amples of phase transitions with singularities appearing 
at various values of A. For example, in a superconduct- 
ing system with the charge transfer quantized in pairs of 
electrons, a singularity would appear at A = 7r/2, (see 
example V.A in Ref.©. 
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